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Continuous wave optically and electrically detected magnetic resonance spectroscopy (cwODMR/cwEDMR) 
allow the investigation of paramagnetic states involved in spin-dependent transitions, like recombination and 
transport. Although experimentally similar to conventional electron spin resonance (ESR), there exist limita- 
tions when applying models originally developed for ESR to observables (luminescence and electric current) of 
cwODMR and cwEDMR. Here we present closed-form solutions for the modulation frequency dependence of 
cwODMR and cwEDMR based on an intermediate pair recombination model and discuss ambiguities which 
arise when attempting to distinguish the dominant spin-dependent processes underlying experimental data. 
These include: 1) a large number of quantitatively different models cannot be differentiated, 2) signs of signals 
are determined not only by recombination, but also by other processes like dissociation, intersystem-crossing, 
pair generation, and even experimental parameter such as, modulation frequency, microwave power, and tem- 
perature, 3) radiative and non-radiative recombination cannot be distinguished due to the observed signs of 
cwODMR and cwEDMR experiments. 
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I. INTRODUCTION 

Electron spin resonance (ESR) is a useful tool for the inves- 
tigation of microscopic properties of paramagnetic states in a 
wide variety of materials. In conventional ESR experiments, 
the total polarization of the investigated spin ensemble is ob- 
served by the measurement of microwave absorption. In some 
materials, there are other observables which can be used to 
detect electron spin states. For instance, when electron spins 
control electronic transitions such as transport or recombina- 
tion, macroscopic materials properties such as photolumines- 
cence, electroluminescence or conductivity can change under 
spin resonance. Fig. Q] depicts a conceptual process of spin- 
dependent recombination 1,2 which can be detected by ODMR 
and EDMR. The advantage of these electrically detected mag- 
netic resonance (EDMR) and optically detected magnetic res- 
onance (ODMR) spectroscopies is that they are significantly 
more sensitive than conventional ESR (spin polarization is 
usually low), and provide direct insight regarding how param- 
agnetic states in semiconductors affect some of the technolog- 
ically most widely used electrical and optical materials prop- 
erties. ODMR has been used in a wide range of research ar- 
eas since its first invention 3 ' 4 . ODMR and EDMR are about 
8 to 9 orders more sensitive than ESR, they both are proven 
to have single spin sensitivity 2 ' 5 " 8 , and they both can directly 
link a paramagnetic center to a specific luminescence cen- 
ter—"^. Thanks to these advantages, ODMR can be used to 
deconvolute unresolved, overlapping luminescence bands in 
semiconductors 10 . EDMR provides information about elec- 
tronically active paramagnetic centers in a similar way, again 
with higher sensitivity than ESR^iii. In the early stage (until 
about the 1980's), ODMR was mainly conducted on inorganic 
semiconductors to identify paramagnetic recombination cen- 
ters and to investigate their spin-dependent processes 5 -*^. It 
played an important role in investigating spin-dependent pro- 
cesses especially in amorphous silicon (a-Si) and revealed a 
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Figure 1 . Spin-dependent recombination via localized paramagnetic 
bandgap states. Excess charge carriers, electrons and holes can re- 
combine via a localized paramagnetic state which acts as a recombi- 
nation center. If a conduction electron and a unpaired electron at a 
paramagnetic recombination center form a spin singlet pair, the con- 
duction electron can be captured by the recombination center. An 
electron at the recombination center can eventually recombine with 
a hole and create a photon. When they form a spin triplet pair, the 
capture probability of the conduction electron by the recombination 
center is low and excess carriers contribute to photocurrent. Because 
recombination process is dependent on mutual spin orientation, this 
recombination rate can be altered by ESR when they are weakly cou- 
pled. Thus ESR can alter recombination rate which results in change 
in photoluminescence and photocurrent, and they can be detected op- 
tically and electrically, respectively. 



variety of defect states which influence recombination in a- 
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Continuous wave ODMR and EDMR (cwODMR and 
cwEDMR) have been used in a wide range of research 
fields: they have been used to investigate spin-dependent 
transitions involving phosphorous donors in crystalline sili- 
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co n 11 ' 17 , trapping centers and their recombination dynamics 
in nanocrystals 8 ' 1 , transport and recombination in micro- 
crystalline hydrogenated silicon 21 , GaN 22 ' 23 , and SiC 2 -, and 
spin-dependent recombination in nitrogen vacancy centers in 
diamon d 25 " 27 . Because cwODMR and cwEDMR can be used 
to distinguish overlapping recombination bands and their dy- 
namics in disordered materials, they have also been used 
to investigate (usually amorphous) organic semiconductors: 
cwODMR and cwEDMR have provided information about 
spin-pairs dominating electronic processes and their transi- 
tions in conducting polymers 28-39 , small molecules 4 ^ 2 -, and 
polymer or small molecule/fullerene blends 43 - 44 . The effect 
of isotopic modification on magnetic field effects in organic 
semiconductors also has been observed by ODMR 45 , and the 
intersystem-crossing time has been extracted from the modu- 
lation frequency dependence 46 . 

Experimentally, cwODMR and cwEDMR are similar to 
conventional ESR except that luminescence intensity and 
electric current are picked up instead of the microwave ab- 
sorption. Two magnetic fields, a static field Bq and oscil- 
lating field B\, are applied to a sample with Bq _L B\. The 
frequency of the sinusoidal B\ field is matched with the Lar- 
mor frequency of the paramagnetic center to satisfy the res- 
onance condition. As for most ESR spectrometers, X-band 
(« 9.7 GHz) is used, a frequency in the microwave (MW) 
range. In the case of cwODMR, to allow for optical detec- 
tion, optical or electrical excitation of electronic states is nec- 
essary. Depending on the excitation method, photolumines- 
cence detected magnetic resonance (PLDMR) or electrolumi- 
nescence detected magnetic resonance (ELDMR) can be per- 
formed. In the case of PLDMR, constant optical excitation is 
applied using, for example, a Laser, and the resulting photo- 
luminescence (PL) is detected. To increase the signal to noise 
ratio, lock-in detection is oftentimes employed. Two differ- 
ent modulation methods can be used. One method involves 
modulation of the static magnetic field, Bq, as used for con- 
ventional cwESR, the other approach is based on the modu- 
lation of the MW amplitude. Experimentally, Bo modulation 
has been found to give weaker signals than MW amplitude 
modulation 5 . Square modulation of the microwaves at a fixed 
reference frequency is generally used. The PL intensity re- 
flecting the varying MW amplitude is then fed into a lock-in 
amplifier, and both in-phase and out-of-phase signals are ob- 
tained. In some studies found in the literatur e 2 1 i 3Q i 38 - 4 1 1 47 ? the 
out-of-phase signal is ignored, however, doing so can result in 
the loss of important information, as will be explained later. 

When the optical excitation is also modulated, a double 
modulated PLDMR (DMPLDMR) becomes possible^. An 
experimental setup for a MW modulated ODMR experiment 
is shown in Fig. [2] For EDMR, the optical detection is 
replaced by a current measurement. The metallic contacts 
needed for this, require a design that prevents the distortion 
of the MW field. 

For both cwEDMR and cwODMR, the responses of the 
observables to the induced magnetic resonances are de- 
termined by the underlying electronic processes. The 
time scales on which these processes occur depend on 
various experimental parameters, such as excitation den- 
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Figure 2. (Color online) Sketch of a cwODMR setup. The basic prin- 
ciple of cwODMR is the same as that of conventional ESR. Square 
microwave modulation can be used instead of Bo ue ld modulation 
and a lock-in amplifier is employed to increase the signal-to-noise 
ratio. 



sity 6,7,i3,i8,4i,48-50 (or an m j e ction current for EDMR 30 i 41 i 51 ), 
temperatur e 6 ! 1Q ' 30 i 41 , and MW power (equivalently B\ field 
sti-ength) 6 ! 7 ! 12 ! 17 " 20 ! 26 ! 40 ! 43 ! 46 ! 52 -.^ 5 .. The dependencies of 
cwODMR and cwEDMR signals on these parameters can al- 
low us to distinguish overlapping transitions and to understand 
their dynamics. For cwODMR, spectral information also can 
provide additional information for distinguishing overlapping 
luminescence bands 6 ,10 ' 26 . 

Another experimental parameter that can influence the ob- 
served cwODMR and cwEDMR signals is the modulation fre- 
quency, as the lock-in detected signals depend on the transient 
responses to the modulated M W 6 ' 10 ' 13 ' 20 ' 56 . Although its im- 
portance has been sometime discussed in conventional ESR 
studies 5 - 7 -^, modulation frequency effects on cwODMR and 
cwEDMR have often been ignored in the literature, and, as a 
result, studies often reported results obtained using only one 
(or a small number) of modulation frequencies (usually the 
one which maximized the obtained signal). One can, how- 
ever, find a number of reports showing modulation frequency 
dependencies. Different signals at different modulation fre- 
quencies were reported for the first time by Biegelsen et al 53 . 
Other investigators have noticed that modulation frequency 
effects play an important role in the observed signal, which 
can change drastically as a function of the modulation fre- 
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Qualitative reports of modulation frequency de- 
pendencies can be found in the early ODMR and EDMR liter- 
ature 5,10 ' 14 which were sometimes used to identify the overlap 
of separate spin-dependent signals 7 . Even so, very little sys- 
tematic research into modulation frequency effects was under- 
taken before the late 1990's, when research into this question 
became more common 1 9 1 20 1 2e 1 3 1 ^37—39,46,49- ^5^ 

A number of researchers have attempted to understand 
modulation frequency effects by developing rate models. 
Dunstan and Davies were the first to develop solutions for 
ODMR transients^ 2 -. Next, Street and Depinna et al. devel- 
oped rate models and found transient solutions&ii Lenahan 
et al. explained their observed modulation frequency depen- 
dence using a simple rate model described by only one time 
constant 15 . A number of studies based on the steady-state 
solutions of such rate models have been reporte d 1 7 ' 1 W 2 i 54 . 
However, to understand the modulation frequency effects the 
exact solutions for the frequency dependence are necessary. 
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There has been a number of efforts to find the solutions for 
modulation frequency dependence 8 ' 20 ' 3 1 i 37 - 39 - 43 i 46 - 50 . How- 
ever, no closed form analytical solutions have been reported, 
and important aspects of modulation frequency effects remain 
not well understood. This has led to a number of debates 
regarding the underlying physical mechanisms of cwODMR 
and cwEDMR signals, because modulation frequency depen- 
dencies observed by different groups on similar systems have 
sometime led to completely different spin-dependent transi- 
tion models. For example, the source of EDMR and ODMR 
signals seen in organic semiconductors has been attributed to 
both a spin-dependent polaron pair mode l 39 ' 46,59 and a triplet 
exciton-polaron quenching mode l 37,38 ' 56 . 

Lock-in detected cwODMR and cwEDMR signals can be 
either positive or negative depending on the shapes of tran- 
sient responses 6 , 1 °' 1 3 i 20 i 56 . A variety of spin-dependent mod- 
els have been developed based on the observed signs of 
cwODMR and cwEDMR signals as well as experimental pa- 
rameters, like pair generation rates, temperature, MW power, 
and modulation frequency. Examples for such studies exist for 
a-si 2 ' 6 - 7 ' 10 ' 12 - 14 ' 16 ' 53 ' 60 , InP nanocrystals^ ., II-IV semiconduc- 
tor a 19 ' 61 , PbL films 48 , nanoparticles 18 , and organic semicon- 
ductor a 30 ' 3 1 ■ 34 - 35 i 37 - 39 1 42 - 44 - 49 i 50 . For instance, it has been gen- 
erally accepted that radiative and non-radiative recombination 
results in enhancement and quenching of cwODMR signal re- 
spectivel y 2 ' 7 ' 1 8,48 i 62 , and all recombination processes and all 
de-trapping processes result in quenching and enhancement 
of cwEDMR signals, respectively 2 ^. The qualitative expla- 
nation for signs of cwODMR signals is as following: spin 
resonance induces mixing between triplet and singlet pairs, 
and because initial states are generally dominated by triplet 
pairs due to the fast recombination of singlet pairs, the num- 
ber of singlet pairs is increased at resonance. Thus, the over- 
all transition rate increases 2 . Some studies even concluded 
that a certain channel is radiative or non-radiative, based on 
the sign of the ODMR signal 6 ' 10 ] 13 ' 60 . The idea here is that 
when a non-radiative recombination process is enhanced un- 
der spin-resonance, the competing optically detected radiative 
channels must be quenched. 

The above examples show how critical it is to understand 
how MW modulation affects the observed cwODMR and 
cwEDMR signals. In this report, we employ the widely 
accepted spin-dependent transition model based on weakly 
coupled electron-hole pairs 63 , and find its closed-form an- 
alytical solutions. We then use this solution to explain 
how a broad range of electronic transitions, including re- 
combination, dissociation, intersystem-crossing, pair gener- 
ation, and spin-flips can affect the cwODMR and cwEDMR 
signals. We show how serious ambiguities related to the 
modulation frequency dependencies can arise, which make 
it difficult to determine the fundamental physical processes 
responsible for the observed cwEDMR or cwODMR fre- 
quency dependence. For example, extensive ODMR stud- 
ies have been conducted on organic semiconductors to de- 
termine their dominant recombination processes. A variety 
of models have been suggested based on the observations 
of the signs of cwODMR and cwEDMR such as the singlet 
exciton-quenching mode l 35 ' 37,49 , the triplet-triplet annihila- 



tion model , the polaron-to-bipolaron decay • , and the po- 
laron pair recombination 39 . We show that in many cases, the 
modulation frequency dependence cannot be used for such as- 
signments, since the sign of these signals can be negative or 
positive for both, radiative or non-radiative processes. 

II. MODELS FOR THE DESCRIPTION OF 
SPIN-DEPENDENT TRANSITION RATES 

The first quantitative model explaining spin-dependent re- 
combination was suggested by Lepin 1 who described a ther- 
mal polarization model which predicted a relative change in 
photoconductivity of less than 1 -6 at 300 K for X-band ESR. 
Microwave frequency and temperature dependencies were 
also predicted. However, it turned out that this model could 
neither explain the signal intensity of more than 10~ 3 that was 
observed in undoped a-Si:H at R.T. 1 , and the very weak de- 
pendencies on microwave frequency 64 and temperatur e) 16 ' 65 . 
These problems were soon resolved by another model devel- 
oped by Kaplan, Solomon, and Mott (KSM model) 63 . In the 
KSM model, intermediate pair states exist prior to a spin- 
dependent transition and the spin pair states may recombine 
or dissociate. In addition, it is assumed that spin pairs in the 
triplet state can be annihilated only when one of pair partners 
is flipped by the spin-lattice relaxation process or the induced 
ESR, pairs dissociate otherwise. Thus, the recombination of 
triplet pairs happens only when they experience a transition to 
the singlet state. 

In the past decades, a number of refinements were intro- 
duced to the KSM model, in which spin-spin interactions such 
as exchange and dipolar interactions exist within the pair, and 
spin-orbit coupling that is weak but not negligible is permitted 
such that weak triplet transitions become possible 66 . Because 
the intermediate pairs, consisting of two spins with s=l/2, can 
experience spin-spin interactions, the pair eigenbasis consists 
in general of two parallel states (|7+) and |7L)) and two mixed 
states (| 2) and |3}) which change continuously from | t|) and 
| |t) to |S) and | To) respectively as the spin-spin interactions 
increase. ESR can induce transitions between the eigenstates 
of weakly coupled pairs such as \T+) O 1 tl), \T-) «-> 1 tl)> 
\T+) *± | tt)> and \ T -) I IT)- Thus, when the spin-spin 
interaction is weak, there can appear transitions among all 
four eigenstates and the transition probabilities are functions 
of the spin-spin interaction strength. Note that transitions of 
I It) I tl) are ESR forbidden but possible due to T\ relax- 
ation, and |2) -H> |3) transitions are possible via mixed relax- 
ation processes. To understand the change of spin pair densi- 
ties by ESR induced transitions, a mathematical approach will 
be given. Boehme and Lips have found the effective changes 
of spin densities by solving Louville equations describing the 
propagation of a spin ensemble during an ESR excitation 66 . 
The corresponding Hamiltonian is 

H = HsgaSa + llBgbSb - JS h -S b -D d [3S z a S z b - S b ■ S b ] + Hi 

(1) 

where the first two terms correspond to the Zeeman terms of 
two pair partners, the third and fourth represent the exchange 
and dipolar couplings, respectively, and the last term is the 
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alternating magnetic field. To describe the weakly coupled 
spin pair, the exchange and dipolar coupling constant, J and 
D d respectively, are assumed to be smaller than the Larmor 
separation. The solutions (density matrix elements) for the 
corresponding Liouville equation can be found elsewhere 66 . 
The density changes of each spin state are then given by 66 , 



energy 



Pl4 (t) = p!> u (t), 

p 2 .3(T)=p 2 3 A V (T)±p 2 3^A W (T) 



(2) 



where indices 1 and 4 represent the states \T+) and |7L) re- 
spectively, pf is the initial density, J and D are the exchange 
and dipolar coupling constants respectively, (Oa represents the 
half of the frequency separation between the states |2) and |3). 
A"(t), A v (t), and A"'(t) represent the ESR duration time (t) 
dependencies. When the Larmor separation (which is the dif- 
ference of the two Larmor frequencies within a pair) is larger 
than the applied B\ field strength, only one pair partner can be 
flipped. In this case the T-dependencies become, 



A"(T) 
A"(t) 



a 2 

l-A(T), 

A w (t)=0 



= A(t), 



(3) 



where H = 2nfn a hi represents the Rabi frequency of the 
flipped pair partner. Therefore, the density changes of each 
eigenstates become 



Pi, 4 (t) = p 1 ° 4 (1-A(t)), 
P2,3(t) = P 2 ° 3 A(t). 



(4) 



Because either one of the states 2 or 3 is always involved in a 
possible transition among four eigenstates, any transition will 
cause a decrease or increase of p 2 or p 3 . Density changes in 
state 2 and 3 are equivalent to density changes of singlet and 
triplet pair states. Therefore we don't need to deal with four 
state problems, instead two pair densities of singlet and triplet 
pairs are enough to describe recombination processes as long 
as any coherent spin motion is not of interest. Note that this 
is a valid statement because modulation frequency is typically 
not faster than the time scale of coherent spin motion so that 
all coherent phenomena will be averaged out. This is also 
the reason why all off-diagonal elements p !; - for i ^ j of the 
Louville density matrix can be neglected. Therefore, only the 
singlet and triplet pair densities, n s and n,, will be considered 
in the following section. 

An illustration of the resulting spin pair rate model is given 
in Fig. [3] Prior to a spin pair transition to a singlet state, it 
is in the intermediate pair state. This pair is created with a 
certain rate, G s for a singlet pair and G t for a triplet pair. 
If this process is due to optical generation of electron-hole 
pairs and spin-orbit coupling is infinitely small, G t can be 
considered to be infinitely small. In the other case, if pair 
generation is achieved due to electrical injection of an elec- 
tron and hole, G t /G s becomes three, because a pair will be 
created with a random spin configuration. The pair can re- 
combine to a singlet ground state with a recombination rate, 
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Figure 3. (Color online) The intermediate pair recombination model 
(KSM) as relevant for cwODMR and cwEDMR. Triplet and singlet 
pairs are formed with two constant generation rates G t and G s re- 
spectively. Those pairs can dissociate into free charge carrier states 
with certain probabilities d t and d s (dissociation rates) or can recom- 
bine to singlet ground state with recombination rates r, and r s . A spin 
mixing process can be introduced by ESR externally and this rate is 
described by a. Another spin mixing process, intersystem-crossing 
is described by kjsc- Note that n t and n s represent triplet and sin- 
glet pair densities, respectively. They do no necessarily correspond 
to eigenstate densities. 



r s for a singlet pair and r, for a triplet pair. This pair may 
dissociate into two free charge carriers without recombina- 
tion. This happens at a dissociation rate, d s for a singlet pair 
and d t for a triplet pair. Before a pair recombines or dissoci- 
ates, it can change its spin configuration from singlet to triplet 
or vice versa. This transition is possible via two spin mix- 
ing processes. One is intersystem-crossing, which is equiv- 
alent to a longitudinal spin relaxation process which can be 
defined as a "radiationless transition between two electronic 
states having different spin multiplicities" 67 . Among many 
processes, the spin-lattice relaxation is one of them which 
can cause the intersystem-crossing. The intersystem-crossing 
rate is described by k/$c- The other process is ESR induced 
spin-mixing as can be seen from eqs. (0) and (|4). This ESR- 
induced transition rate is given by a which is proportional to 
the microwave power (c< B\) and dependent on the spin-spin 
interaction controlled oscillator strength of the pair 68 . 

In the following section, a large number of quantitative 
models will be tested with analytical solutions for the observ- 
ables of cwEDMR and cwODMR. Using realistic values for 
each transition probability, we consider experimentally rele- 
vant values for the cwODMR experiment. A wide range of 
transition rates have been reported. Examples include PL life- 
times in a-Si which span 11 orders of magnitude from 10~ 9 s 
to 10 2 s^,; bound pair decay (e-h pair dissociation) life times 
of 5 x 10~ 5 s in polymer-fullerene blends 70 ; fluorescence life 
times of 2 x 1 0~ 7 s and phosphorescence life times of 1 0~ 4 s in 
conjugated polymers 71 ; microsecond-millisecond time scales 
of recombination in nano-crystalline Ti0 2 thin films 72 ; radia- 
tive decay rates of 10 6 ~ lO 7 * -1 , non-radiative decay rates of 
10 9 ~ 10 10 i _1 , dissociation rates of 10 7 i _1 in organic semi- 
conductors 73 , and a lower limit for the intersystem-crossing 



5 



time of 10~ 5 s in organic semiconductors 74 . In the following 
work, we vary the electronic transition rates, including recom- 
bination, dissociation, intersystem-crossing, and flip-flop, in 
the range between 10~ 4 and 10 9 i _1 to cover as wide a range 
of experimentally observed parameters as possible. 



I 

dn s _ 

dn t 
dt 

where F is the Fermi-Dirac distribution function, F = (1 + 

cir) , which approaches one at low temperature and 1/2 at 
high temperature and used to consider thermalization 18,20 - 46 . 
AE has the order of Zeeman splitting. We chose F to be 
0.25 in all numerical calculations to describe the two-level 
spin system which represents neither a complete thermaliza- 
tion nor a complete non-thermalization. It should be noted 
that a is turned on and off for each half cycle because of the 
square modulated microwave with frequency of 1/T. C s and 
C t are singlet and triplet pair annihilation rate coefficients re- 
spectively. They consist of recombination and dissociation 
rate coefficients, C Sj t = r Sjt + d s<t . Some aspects with regard 
to radiative and non-radiative recombination rate coefficients 
should be mentioned: For radiative recombination, the spa- 
tial correlation between the electron and the hole influence 
the transition probability, so r t and r s depend on the separa- 
tion between an electron and hole 21 ' 7 -. Therefore, because 
the higher generation rate results in less separation, the radia- 
tive recombination probability is also a function of the gen- 
eration rate. However, this effect will not be considered in 
this study, as we assume that the average separation is larger 
than the localization radii of electrons and holes. Note that 
this transition corresponds to the radiative tunneling in hy- 
drogenated amorphous silicon 75 . Non-radiative recombina- 
tion includes all recombination processes which are not me- 
diated by emission of photons, but phonons and hot carriers: 
phonon emission, Auger processes, surface and interface re- 
combination, and recombination through defect states 76 . Non- 
radiative processes quench radiation efficiency in both organic 
semiconductors 77 and inorganic semiconductors 76 . As treated 
by List et al^ and Dyakonov et al4i, we consider both ra- 
diative and non-radiative recombination processes, and thus 



A. Rate equations 

CwODMR is fundamentally similar to conventional ESR 
spectroscopy - the one major modification is that the observ- 
able of ODMR is not the magnetization but the change in the 
number of photons induced by ESR. Generally, lock-in de- 
tected modulation of the Bo or the B\ field is used to enhance 
the resulting ODMR signal. For B\ field modulation, square 
modulated microwaves are continuously applied, and the re- 
sponse to this excitation contains various harmonic frequency 
components. In the following we will focus on this kind of 
experiment. 

Based on the rate model described in Section [TT1 two cou- 
pled rate equations for the singlet and triplet pair densities can 
be written, 



(5) 



(6) 

I 

C s = (r s + r s>nr + d s ) and Q = (r t + r t , nl - + d t ) where the sub- 
script nr indicates non-radiative recombination. 

Given the above definitions, the luminescent intensity and 
electric conductivity become 

/ oe r s ti s + r t n t , (7) 

and 

<7 « d s n s + d t n t , (8) 

respectively. We note that electric conductivity is also de- 
termined by the carrier life time and mobility but ignored 
them because they are merely multiplied to the total disso- 
ciation rate (right term in eq. (O) 74 so that will not affect 
the time dependence nor modulation frequency dependence. 
Non-radiative recombination behaves as a pair annihilation 
process as other radiative recombination and dissociation, but 
it does not appear as proportionality constants in eq. (jTJ and 
©. In the following sections, only radiative recombination 
will be considered (r s nr , r t nr = 0) for simplicity and the con- 
tributions of non-radiative recombination will be discussed in 
Section rvTIl It shall be noted that there are many more compli- 
cated scenario for ODMR detected spin-dependent transitions 
conceivable, including ODMR signals due to non-radiative 
spin-dependent transitions which compete with non-spin de- 
pendent radiative processes. The stochastical description of 
these processes with rate equations is more complex but leaves 
the conclusions made in the following for directly detected ra- 
diative spin-dependent processes unchanged. 

Rate equations similar to eq. ((5]) and (|6]l can be found 
throughout the literature. However, usually only steady 



-C s « s + a(n t -« s ) -kisc(n s ~Fn s ) + fcisc(«t - (1 -F)n x ), 



-C t n t + a(n s -n t ) -k isc (n t -(l -F)n t ) +£isc(«s -Fn s ), 
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state solutions were found for the consideration of cwODMR 
and cwEDMR experiment s - 1 7 i 4 V 8 . In some cases, only the 
time dependence was considered 12 ' 13 ' 43 . Modulation fre- 
quency dependence solutions have also been reported, but 
there have been no reports of closed-form analytical solutions. 
Some solutions reported in the literature were obtained from 
a simplified rate model 1 5 i 31 ' 50 , some solutions were based 
on the steady state 37 ' 38 , some solutions based on the rate 
model reported here were solely reported as numerical solu- 
tion a 8 ' 20,39 ' 46,59 , or the described observable was not the num- 
ber of photons or electrons but total spin densities 31 ' 39 ' 46 . One 
solution given by Hiromitsu et al. was based on an assumed 
steady state for the half cycle where the MW is off 4 - 3 -. 

The rate equations corresponding to eq. (0 and © are 
solved for the two separated time regions where the pulse is 
on and off, and the closed-form solutions can be explicitly ex- 
pressed as: 

n.1 (0 = A n e- mi " + A 21 e-"*' + n° sl , (9) 



tj (C s + W\j) (C t + Wzj) - W\ jW 2 j ' 

wn = a + fcisc(l -F), w 2 i = a + k lsc -F, 

wi2=kisc{l-F),w 2 2=k lsc -F, (17) 

where j=l or 2. It should be noted that the exponents, m\j 
and ni2j, are independent on either the generation rates or the 
modulation frequency. It can be easily seen that m.2j is de- 
cided by the fastest rate coefficient, but it is difficult to pre- 
dict mi/. However, it is clear that m 2 j is always larger than 
m\j. Two constant terms, n®: and are the steady-state so- 
lutions which the system assumes for very low modulation 
frequenc y 17 ' 39 ' 46 ' 49 ' 78 . It should also be noted that the singlet 
and triplet pair populations will approach values at the end of 
each half cycle which are at the same time the initial values of 
the following half cycle. Therefore, the frequency dependence 
can be explained in terms of the differences between the pop- 
ulations at the end of each half cycl o 39 ' 46 , n s \(T /2) — n&iT) 
and n t i(r/2) — na(T). However, lock-in detected signals are 
not simply decided by these quantities. The observables are 
not the population changes, but the changes in the number of 
photons, which incoprorates both the population change and 
the recombination probability. 



n a {t)=Bue- mut +B 2 ie- m ^ +n° tl , (10) 



n S 2(0 =A 12 e~ mu( - t ~ 1 z ) +A 22 e~ m22(r ~T ) +n% (11) 



n t2 (t) = B l2 e- mi2{ '-^ +B 2 2e-' n22{ '-^ +n%, (12) 

where n s \ and n x \ are the singlet and triplet populations when 
the MW pulse is on, and n S 2 and n t 2 are the singlet and triplet 
populations when the MW pulse is off. Those solutions con- 
sist of double exponential functions as is often found in the 
literatures regarding pulsed experiment s 66 ' 74 ' 79 ' 80 . 

The introduced constants in the above solutions are sum- 
marized below, 



(13) 



(14) 



Because the spin populations assume the steady state only 
as the modulation frequency / — » 0, the time dependent so- 
lutions must be solved to explain the transient behavior at ar- 
bitrary modulation frequencies. To find the exact solution, 
the expressions for the eight unknown coefficients Ay and Z? i; 
(/, y' = l orl) in eq. (0, ( [Tol l, ( fTTT i. and (fT2l must be derived 
by application of eight boundary conditions. The boundary 
conditions used as well as the subsequent derivation of the 
analytic form of the coefficients are given in Appendix [Al 

Equations (IA14b . (IA15t . (!A161 l. (IA17b represent exact and 
general analytical solutions for the singlet and triplet den- 
sity functions during a cwODMR modulation cycle. We are 
thus in a position to determine the temporal evolution of the 
cwODMR observable. 



HI. TRANSIENT BEHAVIOR OF CWODMR 

The observable in cwODMR is the emission rate of pho- 
tons, and, as described in eq. (jTJ, the time dependence can 
be obtained by adding the contribution from the singlet and 
triplet pair populations multiplied by the singlet and triplet re- 
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Figure 4. (Color online) A time transient calculated from a numerical 
model described by a combination of parameters as r s = 10 4 
r t = 10° s~\ d s = 10 2 s~\ d t = 10 6 a -1 , £ ISC = 1(T 2 s' 1 , 
10 5 f _1 , F = 0.25, G s = 10 23 f _1 , and G, = 10 20 



-l 



s 

S~ 1 , 0! = 

The dash- 



dotted curve shows the overall response obtained from eq. j I St and 
l |19l >. The blue solid and red dashed curves are the in-phase and the 
out-of-phase components described by / s i sin( and / c i cos(?j£?), 
respectively. See detail in text. 

combination rate coefficients respectively. Thus, 

h = (r s An+r t B n )e- mn < 



(18) 



h = {rAn + nBn)e- mil{ '-^ 
+ (r s ,A22 + r t B 22 )e- m ^'-^ 



's«s2 



r t n° a 



(19) 



where, I\ and l 2 are the photon emission rates due to recom- 
bination of both singlet and triplets pairs when the pulse is 
on and off, respectively. The dash-dotted curve in Fig. |4]is a 
numerical example of the time dependence. Because mi; and 
m 2 j are always positive and m 2 j > mi;, the first and second 
terms in both eq. (fT8l and ( [T9| > determine the slower and faster 
decay, respectively. It is difficult to predict which response 
will show an enhancement or quenching behavior because the 
overall response depends not only on mu and m 2 j but also on 
r s A{j + r t Bjj. Since the coefficients of all exponential terms 
have very complicated dependencies on a variety of parame- 
ters (see eq. (1A14I >. ( |A15t , (IA16I ), and dA17| )), it is clear that 
sign predictions depend on the magnitudes of many parame- 
ters at the same time. Using the above solution, we have been 
able to reproduce a wide variety of cwODMR transients re- 
ported in the literatur o 6 i 10 i 12 i 13 ] 20 . 



IV. MODULATION FREQUENCY DEPENDENCE 

The time dependence solutions, eq. ([TBI and ([T9i l, are the 
collective responses to the modulated B\ field over all fre- 
quency ranges. However, in experimental implementations 
which utilize a lock-in technique, only the component of the 
transient signal which has the same frequency as the refer- 
ence will be obtained. With lock-in quadrature detection, 
both the in- an out-of-phase components are available. While 
the out-of-phase components have often been ignored in the 
literatur o 2 1,30 ' 38 i 4 1 1 47 , we note that the out-of-phase compo- 
nents contain important information, which has been some- 
time addressed in longitudinally detected electron spin reso- 

CT CO 

nanc o J . 

The details of the modulation frequency dependence solu- 
tions are given in Appendix IB1 We are thus able to find an 
analytic expression for the in-phase and out-of-phase com- 
ponents of the transient during a MW modulated cwODMR 
experiment. These are given by 



Vin 



Vb 



1 



— cos(pi) = -I sU 



V ut= ^y-sin(<pi) = i/ c i 



(20) 



(21) 



where Voi is the magnitude of the first harmonic component, 
7 s i and 7 c i are the amplitudes of the first sine and cosine com- 
ponents, and <pi = tan -1 ( ! f L ) (see AppendixlBl. 

Thus the in-phase and out-of-phase cwODMR signals are 
the Fourier coefficients of the lowest frequency sine and co- 
sine terms of the Fourier series solution (eq. (IB1I) ). respec- 
tively. Examples are shown in Fig. [4] to explain the decom- 
posed in-phase and out-of-phase components of the time re- 
sponse. It should be noted that the cwEDMR solutions can 
also be obtained by replacing r s and r t in front of the exponen- 
tial functions in eq. ( fT8l and $1% with d s and d t respectively 
as shown in eq. (O 79 . 

Similarly the solutions for Z?o-field modulated cwODMR 
and cwEDMR can be found in the same way as for microwave 
modulated cwODMR and cwEDMR. While the difference be- 
tween these two modulation techniques is that the spin reso- 
nance is modulated by a square function and a harmonic func- 
tion respectively, the lock-in detected observables are identi- 
cal since the lock-in technique is sensitive to the lowest har- 
monic component only in either case. 



A. At low modulation frequency 

We use the low modulation frequency limit to check the 
solution of our model, by varifying that these solutions can 
explain the cwODMR response. From the solutions above, 
the low frequency behavior is seen to be 
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_ _ (r s + r t )(G t + G s )a + (r t r s + r s w 22 + r t wi 2 )(G t + gs) + rAgt + MtG s _2 
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Figure 5. (Color online) Three different quantitative models result in 
indistinguishable frequency dependencies. Each quantitative model 
is determined by a different set of parameters. Refer to Table|I]for all 
used values. 



V, 



out.lf : 



0. 



(23) 



The out-of-phase component vanishes since the transient re- 
sponse can easily follow the slow modulation. The in-phase 
component shows a typical microwave power dependence: it 
vanishes at small power (when a — > 0) and it becomes sat- 
urated at high power (i.e. it has a non-zero constant value). 
The MW power dependencies of eqs. d22i i and (l23l will be 
explained in the Section [Vl 



identical frequency dependencies of three very different quan- 
titative models. The frequencies at which the in-phase signals 
have their maximum slope and the out-of-phase signals show 
their local maximum values are almost identical, and their 
shapes are also indistinguishable. The patterns shown in Fig. [5] 
represent in fact the most common frequency dependency that 
we have found out by the tested quantitative models. This il- 
lustrates the difficulty in extracting correct values for the cor- 
responding coefficients from a simple frequency dependence 
- one can find a wide range of values which can reproduce 
it. This ambiguity is one of the most significant disadvan- 
tages of cwODMR or cwEDMR. It puts many interpretations 
of cwODMR data reported in the literature in question. 

Of the nearly thousand models we tested, we were able to 
describe them all with only seven frequency dependency pat- 
terns. These are shown in Fig. [6] We find that those pat- 
terns are determined mostly by the recombination rate coef- 
ficients, the microwave power, the spin mixing rates, as well 
as the generation rates. How each parameter influences the 
frequency dependence will be discussed in the following sec- 
tions. The most trivial cases, seen in Fig.[6](a) and (c), will be 
discussed first. 



C. Trivial case (small spin mixing rates) 



B. Ambiguity of cwODMR measurements 

To understand the modulation frequency dependence of 
cwODMR, we inspected a large number of quantitative mod- 
els. There is an extremely large number of possible qualita- 
tive and quantitative relationships betwen the model param- 
eters. To limit the number of cases that we inspected, we 
choose a number of relationships between these parameters. 
We considered that i) the triplet recombination coefficient is 
the smallest one among all the recombination and dissocia- 
tion rate coefficients (r t < r s , d s , c/ t ) (unless otherwise noted), 
and ii) the singlet dissociation rate coefficient is smaller than 
the triplet dissociation rate coefficient (d s < d t ) which means 
that the singlet intermediate state is assumed to be energet- 
ically lower than the triplet intermediate state (unless other- 
wise noted). Under these assumptions, a large number of 
quantitative models were investigated by varying r t , r s , d s , d t , 
k\sc, and a in the range from 10~ 4 to 10 9 s~ l . We investigated 
almost a thousand different variations of the relationship be- 
tween those parameters. 

After looking through these cases, we find that it is al- 
most impossible to distinguish some of the quantitative mod- 
els based on their modulation frequency behaviors. Fig. [5] il- 
lustrates this ambiguity. Figure[5](a), (b), and (c) show nearly 



To understand the behavior of the response to the modula- 
tion frequency, the trivial patterns will be discussed. "Trivial" 
means that the spin mixing rates, both &isc and a are negligi- 
ble when compared to all the other rates. In this case, only the 
spin pair annihilation processes determined by the recombina- 
tion and dissociation rate coefficients become dominant. All 
the patterns in Fig.|5]as well as the patterns in Fig.|6](a) and (c) 
are obtained under the assumption of insignificant spin mix- 
ing rates, %c and a. The pattern in Fig. 0(c) is identical to 
the one in (a), but inverted due to different ratios between G s 
and G t . We found that the sign of the lock-in detected signal 
depends on almost all transition processes as one can deduce 
from Table U 

The most often seen patterns displayed in Fig.|6](a) and (c) 
can be described qualitatively as following: at low frequen- 
cies, the in-phase signal has a constant non-zero value with 
no out-of-phase component. This is because the approach to 
the steady-state takes place on a time scale much faster than 
the modulation period, and the recorded transient response 
looks like the applied microwave pulse train shown in Fig. UJ 
(a). The in-phase and out-of-phase responses are not signif- 
icantly changed until the modulation frequency approaches 
the slowest time constant, mT, , as one can see from the low- 
frequency responses in Fig. [5] For all cases in Fig. [5] and in 
Fig. [71 m\j and m 2j are in the ranges of 10 2 s ~ 10 6 s~ l and 
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Figure 6. (Color online) Seven distinguishable patterns of the mod- 
ulation frequency dependence of cwODMR have been found out of 
almost a thousand quantitative models, (c), (d), and (f) are equiva- 
lent with (a), (b), and (e), respectively, but with opposite signs. The 
parameters used for this data are listed in Table|I] 



10 4 . 



~ 10 6 s , respectively. As the modulation frequency 
approaches m\s, the system begins to lag behind the applied 
MW modulation, and the overall response ceases to resemble 
the simple harmonic function. This results in a decrease of the 
in-phase signal and an increase of the out-of-phase signal as 
seen in Fig.[7](b). At very high frequencies, much faster than 



than the fastest time constant, m 2 J 



lfr 



, both the in- 



15,50 



and out-of-phase components tend to approach zero. This be- 
havior is explained by the exponential decay functions which 
become linear with small arguments and thus, they become 
constants (no change) when the period, T 



D. Recombination, dissociation, and flip-flop 

Because cwODMR measures emission rates of photons, 
which are usually determined by the dominant singlet recom- 
bination rate r s n s , one might expect that r s has a dominant role 
in determining the frequency dependence pattern. In general, 
this is not the case though: other rate coefficients, especially 
spin mixing rates, can dominate the behavior of a cwODMR 
signal. Fig. [8] shows one of the most frequently observed ex- 
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Figure 7. (Color online) Calculated transient behaviors at different 
modulation frequencies. Black dash-dot line is overall response and 
blue solid line and red dashed line are in-phase and out-of phase 
components of it. Parameters are the same with those in Fig. [5] (a). 
The three graphs are normalized by the same scaling factor. Thus the 
relative intensities among three graphs can be compared. 
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Figure 8. (Color online) Role of the singlet recombination rate, r s . 
When r s is small, no significant change in the frequency dependence 
pattern is found when a is increased (from (a) to (c)). But for large 
r s , a pattern change is observed when a is increased (from (b) to 
(d)). All four quantitative models have the same parameters except 



(a)r s = 10 2 J" 1 ,a : 

r s = 10 2 s _1 , a = 10 5 s -1 , and (d) r s = 10' s~\a= 10° s~ l . The 
values for the other parameters used for this data are listed in Table 

m 



1(T 3 s-\ (b) r s = 10 7 s~ 



,(c) 



amples of the frequency dependence patterns influenced by 
both r s and a. 

When a is small, increasing r s has little impact on the ob- 
served frequency dependence (Fig. [8] (a) and (b)). The most 
significant effect is a shift of the frequencies where both the 
in-phase and the out-of-phase components show their maxi- 
mum rate changes. This is due to the increase of the time 



constants, m^, from m\j ~ 10 4 s 1 and m,2j 



10 6 s- 1 to 



m\j ~ 10 6 s 1 and ni2j ~ 10 7 s due to very fast r s . It 
should be noted that d t is 10 6 s~ l in all examples in Fig. [8] 
and r s is 10 7 s~ l in Fig.[8](b) and (d). The frequency depen- 
dence also shows little change when r s remains small and a is 
increased (Fig. [8](c)). This corresponds to Fig. [6] (a) and (c). 
However, when a becomes fast enough to compete with the 



10 



Increasing r 





— I — I — I — I — I — I — I — I — I — I — 

(a) 




1 1 1 1 1 1 1 1 1 1 

(b) 


'E 

3 0. 




urb. unil 

o 










f 


/ l r 10 ' 






c 




c 






k , s r 10 " 




k |K =10 


io 1 io 1 io 3 io s io 7 io 9 


itf* io' io 3 io 5 io 7 io 9 




modulation freq. (Hz) 




modulation freq. (Hz) 


i. units) 


— I — I — I — I — I — I — I — I — I — I — 

(C) r=10 ! 




— I — I — I — I — I — I — I — I — I — I — 

(d) ~ r,=10' 


^N*^ k isc =1 ° 8 


'E 

3 


/ \ k sc =10* 






-Q 

•f: o- 




S n. 








intensi 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — , 


c 

<u 
c 






10 10 10 10 10 
modulation freq. (Hz) 

In-phase 



modulation freq. (Hz) 

Out-of-phase 



Figure 9. (Color online) Role of the intersystem-crossing rate, A^gc- 
At small r s , local extrema appear on both, the in-phase and the out- 
of-phase signal at the high frequency region, when kisc becomes 
large (from (a) to (c)). At large r s , the in-phase signal shows local 
extrema and the out-of-phase signal shows sign change as Aqs C is in- 
creased (from (b) to (d)). All four quantitative models have the same 



(b) r s = 10' ff-\ fc ISC = IO" 2 *" 1 , (c) r s = iu- s ,Kjs C - 
(d) r s = 10 7 s _1 , feisc = 10 8 s — 1 . The other parameters used for this 
data are listed in Table[TJ 
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slower time constant, nf[j, (or even faster than m^), and r s 
is faster than any dissociation rate coefficients, a more com- 
plicated frequency dependence emerges. The in-phase signal 
now has a local extremum. The out-of-phase signal not only 
shows the local extremum (as in the simple pattern) but also 
a zero-crossing point, due to a sign change (Fig.[8](d)). This 
pattern corresponds to Fig. [6] (b) and (d). In this section the 
intersystem-crossing rate, Aqsc, was chosen to be small to in- 
vestigate the influence of a. We note that this pattern also ap- 
pears when kisc becomes large with a small a, as explained 
further in the following section. Note that for cwODMR ex- 
periments this pattern appears only when r s becomes faster 
than any dissociation rate coefficient and a or Ajsc is fast, too. 
It can also be seen for cwEDMR experiments when the dis- 
sociation rate coefficients and a or Aqsc are fast (not shown 
here). We can thus infer that dissociation has a similar effect 
on cwEDMR experiments as recombination on cwODMR ex- 
periments. 



E. The influence of intersystem-crossing on cwODMR 
experiments 

Because the intersystem-crossing rate, k\sc, represents a 
spin mixing process, it acts in a similar way as a, even though 
a is modulated in time. To investigate the influence of fcjsc, CC 
is chosen to be small in this section. When kisc is slow, very 
little change of the frequency dependence as a function of r s is 
seen, similar to the behavior described in the previous section 
(Fig. [9] (a), (b)). In contrast to the case of large a and small 



r s , a major change in the frequency dependence can be seen 
at fast k[$c an d slow r s (Fig.[9](c)). A second local extremum 
appears in the out-of-phase component and a small bump at 
high frequency in the in-phase component. When both &isc 
and r s compete with each other, a new pattern appears (Fig. [9] 
(d)). This pattern is similar to Fig.[8](d) and similar to the pat- 
tern in Fig.|6](e) and (f) when Vi n i f — > at small a (eq. (122b). 
The other important observation is that the sign changes from 
positive for both in-phase and out-of-phase components ((a) 
and (b)) to negative ((c) and (d)). These sign changes due to 
kisc is explained in sec. I VIII 



F. Pair generation 

Due to spin-selection rules, optically generated electron- 
hole pairs (the geminate state) are formed in singlet states and 
remain in this configuration unless strong spin-orbit coupling 
is present 81 . Thus, we can assume G s 3> G t . Figure [6] (a) 
corresponds to this case in which the in-phase and the out-of- 
phase components are always negative and positive, respec- 
tively. This case represents the frequency dependence of pho- 
toluminescence detected ODMR (PLDMR). In contrast to op- 
tical generation, the spin configuration of electron-hole pairs 
formed electrically, i.e. via electrical injection, is determined 
by spin statistics and we can assume 3G S f=s G t . All parameters 
in Fig. |6](a) and (c) are the same except that G s = 10 4 x G t 
in Fig. [6] (a) and 3G S = G t in Fig. [6] (c). We can see from 
these calculations that electroluminescence detected ODMR 
(also called ELDMR) can show the opposite sign compared 
to PLDMR, for very similar underlying physical processes. It 
should be noted that this inversion could be found only for cer- 
tain parameter sets, and this inversion can also happen when 
3G S ^ G t - For example, the sign of the in-phase component 
also becomes positive (not shown here) if every parameter re- 
mains the same except for G s = 10 x G t . Thus, cwODMR can 
result in a positive in-phase and a negative out-of phase signal 
even though G t is smaller than G s . This is because the sign in- 
version is also determined by rate coefficients and not just the 
generation rates. These cases will be discussed in Section [VTil 



V. POWER DEPENDENCE 

The spin flip rate coefficient, a, is proportional to the ap- 
plied microwave power—. Thus we can calculate the power 
dependence of cwODMR signals. Examples are shown in 
Fig. [10] For low modulation frequencies, (see Fig. [10] (a)), 
a simple saturation behavior is predicted by eq. d22l and d23l . 
Note that the out-of-phase component is not always zero, but 
approaches zero at low frequencies, as expected from eq. d23l . 
The saturation characteristics becomes more complicated as 
the modulation frequency increases. At 10 4 Hz, the in-phase 
component shows a local extremum before it returns to a sat- 
uration value (Fig. [10] (b)). Experimentally this behavior has 
been reported recently for low magnetic field cwEDMR on 
crystalline silicon interface defects 17 . At high modulation fre- 
quency, the in-phase component shows the usual saturation 
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Figure 10. (Color online) MW power dependence. All four quantita- 
tive models have the same combination of parameters except for (a) 
/ = 10 3 Hz, (b) / = 10 4 Hz, (c) / = 10 7 Hz. At low modulation fre- 
quencies, typical saturation curves can be found. At high modulation 
frequency, a non-trivial saturation behavior occurs. Refer to Table 
for the values used for the other parameters. 



Figure 1 1 . (Color online) Example of a modulation frequency depen- 
dence function showing a change from non-zero-crossing pattern to 
a zero-crossing pattern. The only difference between the two quan- 
titative models can be found in the triplet recombination rate coef- 
ficients, (a) r t = 10° s , (b) r t = 10 6 s . Values for the other 
parameters are listed in Table|I] 



behavior (even though its saturation occurs at very higher 
power) but the out-of-phase component shows a local ex- 
ternum before it approaches a saturation value. It also has 
a different sign than at lower frequencies (Fig. [10] (c)). This 
shows that one can find opposite signs of in-phase and out-of- 
phase signals at high MW power and high MW modulation 
frequencies. 

VI. SIGNAL SIGN DEPENDENCIES ON THE 
MODULATION FREQUENCY 

Sign changes of cwEDMR and cwODMR signal have 
been found in InP nanoparticles 20 and organic semiconduc- 
tor s 39 ' 46 . The sign change of the cwODMR response in or- 
ganic semiconductor has been attributed to the imbalance be- 
tween changes in the numbers of singlet and triplet pairs when 
the pulse is on and off, which are equivalent to n s \(T/2) — 
n^(T) and n t \ (T /2) — «t2<T) in our model. The zero-crossing 
point of the modulation frequency dependence function has 
also been used to estimate the intersystem-crossing tim o 39 ' 46 . 
According to those reports, the zero-crossing can appear at a 
certain frequency where the increase of the number of singlet 
pairs is matched with the decrease of the number of triplet 
pairs so that the change in the total number of pairs is zero. 
However, we show here that the zero-crossing can be due to 
not only the imbalance of changes between singlet and triplet 
pairs but also to other more complicated relationships betwen 
physical parameters. 

As can be seen in the solutions of the rate equations given 
above, the frequency dependence is not simply obtained from 
n s i(r/2) —n S 2{T) and na(T /2) — na(T), but has a compli- 
cated dependence on various parameters. Among the quan- 
titative models tested here, zero-crossing behavior is rarely 
seen. Fig. QT| shows one example: no zero-crossing is ob- 
served for small r t , but when r t becomes larger and very close 
to r s , zero-crossing is observed (Fig.QT|(a), (b)). It should be 
noted that the origin of this zero-crossing is not obvious be- 
cause of the complexity of the solutions, although we note 



that « s i(r/2) — «s2(T) and «ti(T/2) — ria(T) do not meet 
each other at the zero-crossing point in this case, in contrast 
to the model described elsewhere 39,46 . Thus the imbalance 
between changes in n s and n t cannot be the reason for the ob- 
served zero-crossing. We note that zero-crossing also can ap- 
pear due to an overlap of two different spin-dependent mech- 
anisms whose signs are opposite (e.g. in cwODMR of a ra- 
diative and a non-radiative channel). Note however that all 
zero-crossing effect demonstrated here resulted from a single 
spin-dependent process. The existence of zero-crossing indi- 
cates that one can observe different signs of cwODMR and 
cwEDMR signals from identical samples at different modula- 
tion frequencies. 



VII. THE INTERPRETATION OF CWEDMR AND 
CWODMR SIGNAL SIGNS 

The signs of the cwEDMR and cwODMR signals have 
long been considered important indicators for the nature of 
electronic transitions. For example, it has been generally 
accepted that radiative recombination results in positive in- 
phase ODMR signal s 2 i 7 i 18 . However, the recent observations 
of sign changes 20 ' 39 ' 46 at certain frequencies suggest that signs 
may depend on complicated processes and the interpretation 
based exclusively on the sign of a modulated cwODMR or 
cwEDMR signal is not possible. 

CwEDMR and cwODMR signal signs are determined by 
the transient responses of optical or electrical observables to a 
repeated change between on- and off -resonance, as described 
in Section [IV] Because the time constants and pre-factors of 
the double exponential functions in eq. (|9), ( TTOb . ( fTTT l. and ( TT2b 
are functions of all the transition rate coefficients, there are 
many scenarios which can produce quenching and enhance- 
ment signals for both radiative and non-radiative ODMR sig- 
nals as well as for EDMR signals. Many transitions compete 
with each other. For instance, recombination as well as dis- 
sociation are pair annihilation processes but only recombina- 
tion causes PL while dissociation does not. Thus, when a ra- 
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d/r^ d/r_ 

Figure 12. (Color online) Sign changes due to various rate coef- 
ficients, (a) In-phase intensities of the zero modulation frequency 
component as a function of Gt/G s and d/r s , To distinguish positive 
values and negative values, different color scales are used (positive in 
upper left corner, and negative in lower right corner). The black dot- 
ted line describes the boundary separating positive values and neg- 
ative values, (b) and (c) are two randomly chosen two dimensional 
subsets of the data in (a) representing a generation rate ratio depen- 
dence and dissociation rate ratio dependence. These slices are shown 
as white dashed lines in (a). Intensities in (a), (b), and (c) are nor- 
malized but in the same scale, (d) Changes in the numbers of singlet 
pairs, — «j 2 as a function of the same parameters as in (a), (e) 
Changes in the number of triplets pairs, — n% as a function of the 
same parameters as in (a). Intensities in (d) and (e) are normalized 
but in the same scale. All calculations in this figure are obtained 
from the same condition of r s = 10 4 s — r t = 1 s _1 , k\sc = 1 ■> \ 
a = 1 s-\F = 0.25, G s + G t — 10 16 j -1 . 



diative recombination process is slow and dissociation is fast, 
the resonant response may lead to quenching. This example 
shows that the following qualitative description of the sign of 
cwODMR signals is important. 

The study of the sign change of cwODMR signals as func- 
tions of all individual parameters is beyond the scope of this 
work. Instead, only the low modulation frequency behavior 
will be discussed. 



A. For the case of radiative recombination 

As mentioned in Section lHAl only radiative recombination 
has been considered so far and the solution for the in-phase 
cwODMR signal when radiative recombination is dominant 
at low modulation frequency is given in eq. (1221) . A quantita- 
tive analysis has been done by calculating V; n if while chang- 
ing some parameters, for the example shown in Fig. [T2l we 
assumed that both singlet and triplet dissociation probabilities 
are not distinguishable, two mixing rate coefficients, feisc and 
a, are slower than any other recombination and dissociation, 
and total generation rate, G s + G t is fixed to 10 16 s , r t to 



1 and F to 0.25. Fig. Q~2] (a) shows the zero frequency 
in-phase cwODMR signal, Vi n i f, as a function of the rela- 
tive ratio of the triplet generation rate to the singlet generation 
rate, G t /G s , and the ratio of the dissociation rate coefficient 
to the singlet recombination rate coefficient which is fixed to 
r s = 10 4 s . Color reflects the normalized intensity of Vj n if. 
It should be noted that positive and negative values are inten- 
tionally placed in different scales to make them clearly distin- 
guishable. One can find two noticeable features, (i) The inten- 
sity tends to increase as G t /G s becomes larger and becomes 
negative at low G t /G s , as in Fig.Q~2](b). (ii) The intensity also 
depends on the dissociation rate coefficients: when d is larger 
or smaller than the singlet recombination rate coefficient r s , 
Vi n if becomes very small, and shows an extremum and sign 
change. Fig. Q~2] (a), (b), and (c) show that the signs are posi- 
tive at high triplet generation rates and low dissociation rates 
or, equivalently, high recombination rates. When dissociation 
is not fast, signs are positive as long as triplet generation is not 
slower than singlet generation rate. This means that changing 
pair generation between optical and electrical can induce a 
sign change in cwODMR. This behavior can be more easily 
understood by consideration of competing singlet and triplet 
pairs. In Fig. Q~2] (d) and (e), the differences — n'L and 
n ti — n t2' aie plotted for the same parameters as (a). Note that 
the low-frequency solution for the in-phase cwODMR signal, 
Vin,if> i s proportional to ft ("si — n s2) + r s( n tt — n o)' Both plots 
show different behavior compared to Vi n if but the boundaries 
dividing positive and negative values are very similar. When 
the pair annihilation is dominated only by singlet recombina- 
tion, one can infer that the number of singlet pairs quickly 
decreases in the steady-state off-resonance condition. Thus, 
the steady-state is dominated by triplet pairs. Consequently a 
resonant MW converts triplet pairs to singlet pairs, it increases 
the number of singlet pairs which results in an enhancement 
of cwODMR signal. 

This qualitative pictures applies to the region where — 
n® 2 is positive and — n® 2 is negative, in the upper left regions 
in Fig.Q~2](a), (d), and (e) for example. In contrast, if the triplet 
generation is too low (G t < G s ), (lower-left corners in 
Fig. [12] (a), (d), and (e)), only a small number of triplet pairs 
forms during the off resonance steady-state, and the steady- 
state at off-resonance is dominated by singlet pairs. In this 
case spin resonance induced changes to the number of singlet 
pairs can become negative. 

The statements above are based on the assumption of low 
k\sc and a. When k\sc becomes larger than the other rates, 
sign changes are observed as in Fig. [9] and patterns of V- m jf 
(not shown here) similar to the pattern in Fig. [12] are found, 
although slight shifts of boundaries dividing positive and neg- 
ative are seen. Similar shifts have been found at different F 
and a. These shifts can be explained by an expression, 

Gs _ (r s + d s )(w 2 2 - W 2 l) +W 2 2W U - W 2 jW\2 ^ 
G t Ot + dt)(wil -W12)+W22W11 -W 2 \Wi 2 

which is obtained from eq. (122b by setting V; n if = 0. This for- 
mula explains that the boundary separating the positive and 
negative values in Fig. Q~2] (a) is dependent on all rate coeffi- 
cients. Consequently, cwODMR and cwEDMR signs also de- 
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Figure 13. (Color online) The sign of cwODMR signals can be nega- 
tive when radiative recombination is dominant as in (a), and positive 
when non-radiative recombination is dominant as in (b). In contrast 
the signs of cwEDMR are not different, (c) and (d). Used common 
values for each rate parameters can be found in Table U (a) and (c) 
r s = 10 4 , r. m = 1. (b) and (d) r s = 1, r, m = 10 4 . 



pend on intersy stem-crossing rate kisc, the temperature (note 
that F is a function of temperature), and the MW power a. 
We note again that sign changes can also occur at a certain 
modulation frequency as already explained above. 



idence. 

We conclude that dissociation, recombination, ratio be- 
tween singlet and triplet generation, intersystem-crossing, 
temperature, modulation frequency, MW power, and the na- 
ture of recombination (radiative or non-radiative) can all de- 
termine the sign of cwODMR signals. 



VIII. SUMMARY AND CONCLUSION 

A set of rate equations based on an intermediate pair re- 
combination model are presented and generalized analytical 
solutions have been obtained. These solutions have been used 
to calculate modulation frequency dependencies of cwEDMR 
and cwODMR signals. We have investigated how experi- 
mental parameters affect these modulation frequency depen- 
dencies which revealed that a large number of quantitatively 
different models show non-distinguishable modulation fre- 
quency dependence patterns. This implies that the interpre- 
tation of cwODMR and cwEDMR experiments can be very 
ambiguous. We further showed that the sign of cwODMR 
and cwEDMR signals depend on most of the rate coefficients, 
as well as experimental parameters such as temperature, MW 
power, and modulation frequency. Thus, there are many vari- 
ables which can reverse the sign of cwEDMR and cwODMR 
signals and consequently, conclusions about the radiative or 
non-radiative nature of an observed spin-dependent transition 
based solely on the sign of an observed spin-dependent pro- 
cess or its modulation frequency dependence is not possible. 



B. For the case of non-radiative recombination 

Finally, we want to address the question of whether ra- 
diative and non-radiative recombination results in opposite 
cwODMR signal signs. We have checked a number of quanti- 
tative models and two examples are shown in Fig. [13] In con- 
trast to all other cases discussed above, the non-radiative sin- 
glet recombination coefficients, r s nr is taken into account. In 
Fig-EKa) and (c), r s nr is assumed to be smaller than r s to sim- 
ulate the modulation frequency dependence in which radiative 
recombination is dominant. In Fig. [13] (b) and (d), r s nr is as- 
sumed to be the larger than r s to investigate the non-radiative 
process. It should be mentioned again that r Sjnr contributes to 
the pair annihilation process but it does not contribute to the 
radiative emission rate term as explained in Section lHAl Note 
that Fig.[l3](a) shows one of the modulation frequency depen- 
dence patterns that are discussed above. The in-phase signal 
is negative even though r s is most dominant because G s 3> G t . 

Fig. Qj] (b) shows a zero-crossing behavior, thus, the in- 
phase component can be positive and negative even though 
r s nr is dominant. In contrast to the cwODMR cases, the signs 
of the cwEDMR in-phase signals are positive in both cases 
as shown in (c) and (d). To summarize, our results show that 
cwODMR signals can be negative and positive for both radia- 
tive and non-radiative recombination processes. Any conclu- 
sion about the nature of a spin-dependent recombination pro- 
cess from the sign of an observed cwODMR signal is there- 
fore speculative, and should be confirmed with additional ev- 
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Appendix A: Boundary conditions and exact solutions for the 
pair densities 

Four of the boundary conditions can be easily found from 
the periodicity of the solution: n v i(0) = n^{T), n r i(0) = 
na(T), ) = n«2(|). ™d n,i(|) = n tl {\). From these 
boundary conditions, we obtain 

A\ +A 2 + n° sl = A 3 e ( - mi2772) +A 4 e ( " m22r/2) + n° 2 (Al) 
B l +B 2 + n? 1 =B 3 e { - mi2T ^+B 4 e ( - m ^ T ^+n° 2 (A2) 
Aie (-»mT/2) +A2e (-m 2l T/2) +n i ^a 3 +A 4 + «? 2 ( A 3) 
B ie (- m ii r / 2 ) + B 2e (- m 2i r / 2 )+n° =B 3 +B 4 + n%_ (A4) 

After each half cycle, the number of singlet and triplet pairs 
are decreased or increased. These changes depend on the 
given rate coefficients: the number of singlet or triplet pairs 
are increased by pair generation, decreased by the dissociation 
and recombination processes, or either decreased or increased 
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by spin mixing. From this condition, the other four equations 
can be found as 

T T fi 

n s i(-)-n sl (0) = G s - + I (w 2 \na-{C s + wu)n s \)dt, 

(A5) 



2 h 



T T f T 

n S 2{T)-n sl {—)=G s —+ I {w 2 2n t2 - (C s + w l2 )n s2 )dt, 

(A6) 



na(-) -n tl (0) = G t — + / (w n n t i - (Q + w 2l )n tl )dt, 



2 Jo 



(A7) 



T T f T 

na{T)-na{—) = G t —+ I (w l2 n t2 - (C, + w 22 )n t2 )dt. 

(A8) 

By plugging eqs. (l9ll-(fT2]i into ( IA5b - (IA8b . we obtain 

Al ( e (-mnr/2) _ j j +A2 ( g (-m 21 r/2) _ ^ 

_ W2iBi-(C s + w n )A 1 (_ mil r/ 2 ) _ ^ 



W21B2 - (Cs + wn)A2 (g (-m 21 r/ 2 ) _ ^ (A9) 



OT21 



A 3 ( e (-"n27'/2) _ l) + A 4 ( e (- m 22^/2) _ 1) 

W22fi3 - (C, + wn)A 3 . ( _ mi2 r/2) _ ^ 
mi 2 

W 2 2g4 ~ {Cs + Wl 2 )A 4 , (- m22 r/2) _ 
^22 



(A10) 



Bl ( e (->»nT/2) _ 1 ) +B2 ( e (-'»2ir/2) _ 
wnAi - (C, + w 2 i)fli ,(_ mu r/2) 



wnA 2 - {C, + w 2l )B 2 ( _ m2ir/2) 



m 2 i 



(-m 12 r/2) _ 



:( 4 



1) 

1), (All) 



l)+B 4 {e 



(-m 22 r/2) _ 



^12^3 - (Q + W2 2 )B3 



mi 2 

W12A4 - (Q + W 2 2)B4 



l( e (-wi2r/2) 



'«22 



-m^T /2) 



1) 
1) 



1). (A12) 



W22«? 2 



Note that G s + W2i«,°i - {C s + w n )n° sl = 0, G s 
{C s + wi 2 )n° s2 = 0, G, + wnnjj - (C, + W21K , = 0, G t + 
w i2 n ° s2 — {Ct + w 2 2) n ^2 = are used here, which are obtained 
from eq. (fl~5T > and (TT6b . 

Solving eq. dAU - (IA4t . ( |A9MA12t , and by introducing the 

parameters fly = An" = n° - n° An? = n° - 

n 7* ^ 

n,j, and )fy = e~ m <iT , we realize that = Ayj3y and four 
simplified equations 



1 1 -712 -722 \ 

011 jS 2 l -012712 "022722 

711 721 -1 -1 I 

0n7n fti72i -Pu -P22 I 




( A i\ 

An? 
An? 
V An? / 
(A13) 



are obtained for A 



'./■ 



Equation ( IA13b is a fully determined system of linear equa- 
tions which can be solved. This leads to the solution, 



A 22 = (((fti - J8n) ■ (An? - 7ll An?) - (An? - /3 n An?) ■ (y 21 - y„)) 

■{{P21 -Pu) ■ {PuYuYn-Pn) - {Pu 7i2 -PuYn) ■ {P21Y21 -0n7n)) 

-{{P21 -Pu) ■ (7n7i2 - 1) - (0n7i2 - PuYn) ■ (721 - 7n)) 

■((021 -0n) • (An?-j8u7uAn?) - (An? - j3„An?) • (j8 2l7 2i - PuYu))) 

/{{{P21 - Pu) ■ (7ii722 - 1) - {P11Y22-P22Y22) ■ (721 - 7n)) 

■((021 - J3n) • (J8n7ii7i2 - Pn) - (PuYn- PnYn) ■ (P21Y11 ~PuYu)) 

-((j3 2 i -J3n) • (7ii7i2 - 1) - (0n7i2 - 012712) • (721 - Tii)) 

•((021 -01l) • (01 l7ll 722 -022)- (011722-022722) " (021721 -01l7ll))), 



(A14) 



A12 = ((021 - 0n) ■ (An? - 7ll An?) - (An? - /3 n An?) • ( 721 - y„) 

-((021-011) ' (711722 - 1) - (011722-022722) " (jll ~ Yu)) ' A 2l) 
/((02 -01l) ' (7ll7l2 - 1) - (011712 - 012712) ' (721 - 7ll)), 



(A15) 
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A 21 = ((An? - flu An?) ■ (fcifcl - Pnfli) 

-(J811712-J812712) • (fcmi -J8n7u) -An 

-(j8l 1722 -022722) • (fcmi -J8ll7ll) -^22) 

/((fti-jSiO-Cfemi-jSiiTn)), (A16) 
An =An 1 °-A2i + 7i2-Ai2 + 722-A 2 2. (A17) 
Appendix B: Modulation frequency dependence solutions 



2 i_g-">nr/2 

/o= ™(rjAn + r ( Bn)( ) 

T mn 

2 1 — e -' n 2lT/2 

+ ™M 2 i + r f B 21 )( ) 

7 W21 

2 1 — e - ,n n T / 2 

+-{r s A u + r,Bn){ ) 

1 m\2 

2 1 — e~ m2lT l 2 
+ = M22 + r,B 22 ) ( ) 

7 OT22 



+0(«ri+«!fl)+M'«fl+*fl) 



1 „0 



(B6) 



To find the in-phase and out-of-phase components at a given 
modulation frequency, it is better to find the Fourier series of 
eq. ( fT~8b and ( fT9b , and the frequency responses will be decided 
from the Fourier coefficients according to the definition of the 
Fourier series, 



The Fourier series in eq. (IB II ) can be simplified by introducing 
Vb = Vlc 2 +h 2 and <p = tan -1 (£■) as below, 

WO^f + X>sin(2Z7r/f + <p), (B7) 
1 1=1 



/f.v(0 - y + E(^' cos (— +^sin(— r)), 



/, = - / /(f)cos(- f)*, 
7 Jo 7 



**= ~ / I(t)sm(—t)dt. 
7 Jo 7 



(Bl) 



(B2) 



(B3) 



Then the obtained two coefficients as well as the zero fre- 
quency component are: 

, 2m u , w l -e- m i' r / 2 cos(/7r), 

2m 2 i, w l -e- m ^ T / 2 cos(l%), 
+ ^(rsA 2 i+r,B 2l )( 2^,2^2/7-2 

+ _ Ml2+ „ B , 2)( m ; 2 4 4W/r , > 

+ ^ M2;+r , B22)( -('")-f-°:f ), 



m 22 2 + 4/ 2 7r 2 /r 2 



4/tz;, w l-e- m " r / 2 cos(/n:). 

7 * = rf-Mll + ^H 2^/1,2^27x2 

4/7T , , „ .,l-e- m2ir / 2 cos(/7r), 
+ F (rAi + rM m21 2 +W/r 2^ 

w cos(/7r)-e- m i2 r / 2 N 

+ -=2-Ml2 + r,2?i 2 )( 2 1 A ,2^2/T2 

4/71. , „ .,cos(/^)-e- m 22 r / 2 

.0 , .a-Ox/COS^)- 1 ' 



(B4) 



-(r,An? + r,An?)(- 



1% 



where f=l/T is the frequency of the square modulation. A 
Lock-in amplifier multiplies the input signal by its own in- 
ternal reference signals, sin(o)if + Bl) and cos(©£f + Bl), to 
detect in-phase and out-of-phase signals, respectively. Thus, 
the in-phase V,„ and out-of-phase V ou t signals are 

V in = |v L sin(co L r + B L ) 
V l Vo 



I [cos ( (21 %f - co L )t + <p - Bl) 



1=1 



- cos((2l%f + co L )t + <p + 8 L )} , 
Vout = ^V L cos(o) L f + 8 L ) 

+ £ [ sin (( 2/7r / + ®i) f + <p + 0l) 

+ sin((2/7T/- a> L )r + <p - 



(B8) 



(B9) 



where Vl is the amplitude of the reference signals. After these 
signals pass through a low pass filter, only the non-AC signals 
will remain. And the frequency of the internal reference signal 
is fixed such that it has the same phase as the external refer- 
ence signal. Thanks to this condition, col ~ 2%f, the in-phase 
and out-of-phase signals become 



Vb, 



Vbi 



cos (<pi) 



1 



V01 . , , 1 
Vom = — sin((pi) = -7 c i 



(BIO) 



(Bll) 



(B5) 



where Vbi = Vb, I s \ = I s , h i = h, and <pi = <p at / = 1, Bl is 
usually set to zero, and Vl is set to 1 . 



Appendix C: Parameters used for calculation 

The parameters used for all data presented in the figures are 
listed in Table|T] 
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Table I. Parameters used for the calculation of some of the data in this report. All values have a unit of s except for F which is unit-less. 
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